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The physics of an electron solid, held on a cryogenic liquid surface by a pressing electric field, 
is examined in a low-density regime that has not been explored before. We consider the effect of 
the pressing field in distorting the surface at the position of each electron and hence inducing an 
attractive force between the electrons. The system behavior is described in terms of an interplay 
between the repulsive Coulomb interaction and the attractive surface- induced interaction between 
individual electrons. For small densities and large enough pressing fields, we find a parameter 
regime where a square lattice is more favorable than the usual triangular lattice; we map out the 
first-order transition curve separating the two lattice geometries at zero temperature. In addition, 
our description allows an alternate static perspective on the charge-density wave instability of the 
system, corresponding to the formation of multi-electron dimples. 

I. INTRODUCTION 

The crystallization of charged particles due to Coulomb repulsion, first predicted by Wigner 1 , has been under 
discussion in the context of 2D electron systems on cryogenic substrates 2-5 and in semiconductor heterojunctions 14 . 
Wigner crystallization has also been studied in three-dimensional trapped-ion systems 15 . The situation in semicon- 
ductor junctions, unfortunately, has been obscured by the presence of disorder 14 . Electrons held on a liquid helium 
surface by a pressing electric field, on the other hand, provided the first clean realization of a Wigner crystal 3,4 . 

For the case of electrons on a cryogenic liquid substrate, the presence of a distortion-prone surface introduces 
additional physics, fn this paper we investigate a previously-unexplored effect of the surface on the electron crystal, 
in a low-density regime that has not yet been probed experimentally. 

Experiments with two-dimensional electron systems on liquid helium are done with a pressing field E± that holds 
the electrons to the liquid surface. One effect of the pressing field, that has been studied in depth by a number of 
authors 5-8 , concerns a regime of density high enough so that one can treat the surface as a uniformly charged sheet. In 
this regime of large electron density n, as one increases Ej_ to a certain value, the surface breaks up into many-electron 
dimples, the size of the dimples being given by the capillary length of the liquid. The dimples themselves form a 
triangular lattice 5-8 . This situation describes a significant part of the E±-n phase diagram. 

At smaller densities, it seems sensible to consider the effect of the pressing field on individual electrons: each 
electron forms a single-electron dimple. (Since the physics determining the shape of these dimples is the same as 
that for multi-electron dimples, they have the same shape, on different scales.) Single-electron dimples have been 
studied previously, especially in the context of forming a self-trapped polaronic state 9 . However, the formation of 
such dimples for each electron should also cause a 'mattress-effect' attraction between any two electrons, and to the 
best of our knowledge the effects of this surface-induced interaction have not been explored. 

The surface-mediated attractive interaction between electrons introduces new physics in the 2D system. One 
possibility is one or more structural phase transitions. Our analysis shows that, in a low-n, high-i?^ segment of the 
E±-n phase diagram, a square lattice is energetically more favorable than the usual triangular lattice. We thus prove 
the existence of at least one structural transition. In addition, some of the traditionally-known surface effects, such 
as the surface-buckling instability, can be interpreted from a fresh perspective using the idea of competition between 
attractive and repulsive interactions. 

Two-dimensional electrons have been studied on several cryogenic systems 2 , e.g., surfaces of liquid 4 He and liq- 
uid 3 He, at 4 He- 3 He interfaces, the interface between solid and superfluid 4 He, etc. In recent times, experimental 
efforts have concentrated on increasing electron densities, and on using thin liquid films, for example for the purpose 
of observing a solid-to- liquid quantum phase transition in the low-temperature, high-density direction 12,13 . Other 
experiments include investigations of the effects of a magnetic field on excitation and transport properties, and the 
study and control of de-coherence 19 for quantum-computing purposes. 

In section II we derive the form of the attractive interaction between two electrons due to surface deformations 
caused by the pressing field pushing the electrons down. In section III we discuss general properties of systems formed 
by a combination of attractive and repulsive forces, and outline the consequences for the system we are describing. 
The energy calculations for the relative stability of square and triangular lattices are outlined in section IV, and 
section V describes details of the numerics and the resulting phase-diagram. 



2 



II. SURFACE-INDUCED INTERACTION 



We consider a system of N electrons, held at positions on the surface of a thick cryogenic (possibly helium) 
liquid substrate by a pressing electric field of magnitude E± perpendicular to the surface. We will calculate the 
surface-mediated interaction from classical, static considerations. 

We will use u(r) to denote the vertical displacement of the surface at point r, as compared to the undistorted 
(flat) configuration. We proceed to write down the energy as a functional of u(r). There are three contributions: a 
surface tension term describing the energy cost due to surface distortion, a pressing-field term describing the energy 
that the electron gains due to vertical displacement, and a gravity term describing the bulk displacement of helium 
accompanying the surface distortion. 

E[u(v)\ — a J d 2 r [l + (Vu(r)) 2 ] 1/2 + eSx^ufo) + d 2 r[u(v)] 2 . 

Here a and p are respectively the surface tension and density of helium. The pressing electric field E± actually 
contains contributions from both the externally applied field and the field due to the image charge formed by the 
helium dielectric. We will be interested in large applied fields, compared to which the dielectric effect is negligible. 
Expanding the surface-tension term to lowest order in Vu, we get in momentum space: 

4 N A 

2S[«(k)]{'i--™> = ^^VkM-k) + eE^Y,^^ + fEW)' 

k i=l k k 

The k are 2D wave-vectors and A is the area; Air 2 ^ k <-> A J cPk. The form of it(k) is now determined by minimizing 
the energy functional. The result is 



eE ± ^ e~ ik r 



"00 = - -r^>.7^2 = > «?(k). (l) 
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Here Iq = y r a jgp is the capillary length of the liquid substrate; it will be the important length scale in all our 
considerations. Also, 

1 Aa k 2 + l„ 2 

is the Fourier transform of the distortion due to a single electron at r^, i.e., the shape of a single-electron dimple, as 
can be easily verified by minimizing the energy functional for a single-electron system. 

The energy of the system is the minimum of the functional £7[u(k)], and can be now written in terms of the u'*'s: 

4 N N 

k i — l 7—1 i,j k 

which separates into diagonal (i = j) and non-diagonal pieces: 

E(r u ...,r N ) = NE^ + J^V{T i -r j ), (2) 

i<j 



with 

E 2Aa +/ ° > ' V{ri *>> Aa ^ k 2 + L 



The first term of cq (2) is an extensive quantity representing the energy of N independent electrons. The second 
term gives the attractive interaction energy between the electrons mediated by surface deformation. We have thus 
obtained the 'mattress'-effect attractive potential to be 

■rrf \ (effO 2 f ,2, cos(k-r) (e£ ± ) 2 

V ( r ) = - ~02^ / d k u2 , 1 ,12 = 7r—Ko(r/l ). (3) 
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Here Kq is the zeroth-order Bessel function of the second type. It has the asymptotic behavior Kq(x) <~ — log(x) for 
x — > and K (x) <~ (-K/2X) 1 ' 2 e~ x for x — > oo. 

Two comments are in order concerning this derivation. First, the electrons have been treated as point objects, and 
this leads to interactions and dimple shapes that diverge logarithmically at small distances. This divergence is cured 
by the fact that the electron wavefunction has a finite spatial width. In this work we will not attempt a calculation of 
this wavefunction width, but will assume that this spatial extension is much smaller than the lattice spacings which 
we consider. This assumption is particularly reasonable in the low-density region, deep in the solid phase, that we 
focus on. 

Second, we have retained only the (Vu) term in the surface-tension energy. This is equivalent to keeping only the 
pairwise interaction between electrons and dropping four-electron and higher terms. For quantitative calculations, 
eq (3) is used only in a very-low-density regime; effects of four-body or higher terms are expected to be small here. 
Our order-of-magnitude estimates in section III concerning moderate-to-high densities makes use only of qualitative 
ideas from the mattress-effect calculation, because in this regime the effect of higher-order terms is expected to be 
important. 



III. EFFECTS OF ATTRACTIVE INTERACTION: GENERAL REMARKS 

The Wigner lattice on a helium substrate is formed by electrons which interact via a Coulomb repulsion and a 
surface-mediated attraction: 

e 2 (eE,) 2 

V(r) = Koui(r) + y surf (r) = i_^^ ( r // ). (4) 

r Zna 

The effect of the surface term has not been considered in detail before, and we shall proceed to do so in the present 
paper. 

The surface- induced attraction can be tuned by controlling the pressing electric field E± . At short enough distances, 
the 1/r function dominates over the logarithmic surface term, and at large r the Coulomb term again dominates over 
the exponentially decreasing attraction. V SUI { can compete with V^ ou ] only at intermediate distances. As E± is ramped 
up, the distance scale at which V^ ur f first becomes comparable to V con \ is r <~ l - 

In general, when microscopic objects interact via attractive and repulsive potentials of different ranges, several things 
can happen depending on the relative strengths and ranges of the attractive and repulsive forces. First, consider the 
case of repulsive forces alone, or short-range repulsion coupled with longer range attraction. This situation tends to 
create microscopic-ordered phases, such as Wigner crystals, vortex lattices and skyrmion lattices 16 . Second, when the 
attractive force dominates at all distances, the system tends to collapse. One example is what happens at the higher 
critical magnetic field H C 2 of a Type-II superconductor: the interaction between vortices of the mixed phase becomes 
attractive citekramer at H = H C 2, and this leads to collapse of the vortex matter so that the system is filled with the 
normal electrons of the vortex cores, and superconductivity is destroyed. And finally, a combination of short-range 
attraction and long-distance repulsion tends to create intermediate-scale order, or "clustering" . Examples are charge 
density waves in solids and quantum hall systems. 

This approach enables us to view the well-known electrohydrodynamic (surface-buckling) instability of this 
system 5 ~ 8 from a novel perspective. For an electron density much larger than Zq 2 , any electron has a large num- 
ber of electrons within a distance less than the capillary length l Q from itself, so that the attractive force due to 
surface distortion acts between a large number of particles. Thus we have an attractive interaction at intermediate 
distances and a long-distance Coulomb repulsion. Therefore when the attractive interaction is ramped up by increas- 
ing the electric field, one can expect from the preceding general discussion the formation of intermediate-scale clusters, 
which are themselves ordered in a regular pattern. 

This is exactly the phenomenon of formation of many-electron dimples that is observed for high-density electrons 
under large pressing fields. Previously this instability has been studied in terms of the excitation spectrum of a 
surface approximated as being uniformly charged 5-8 . In the traditional analysis, one finds that at a certain pressing 
field, the spectrum goes "soft" at wavenumber k <~ Zjj" 1 , indicating the onset of a charge-density instability of this 
wavenumber. The pressing field at which this spectrum softening first happens is given by El = ^[pga] 1 / 2 - {2nne) 2 , 
or Ej_ :=y [16-7T 2 pg a] 11 4 at low densities. 

In our picture of competition between attractive and repulsive forces, the formation of multi-electron dim- 
ples (intermediate-scale order) would occur when the attractive term V sur f starts to dominate over the repul- 
sive term at intermediate or small distances. This viewpoint allows a simple calculation of the pressing field at 
which the surface-buckling occurs: it is the pressing field for which we have V sur f(r — l ) w V cou \(r — lo), i.e., 

Ej_ w {[if (l)r 2 47r 2 p.gcr} 1 / 4 w 1.1 [l67r 2 p.gcr] 1/4 , within 10% of the traditional result. 
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Next we want to concentrate on a lower- density regime, n > 1/Zq, where we find a less dramatic but nevertheless 
important effect of the competition between attractive and repulsive forces. As the inter-particle distance approaches 
lo, the formation of many-particle dimple becomes less likely, since there is no longer "many" electrons within the 
capillary-length scale, and the uniform-smeared charge approximation becomes untenable. It is not clear whether a 
surface-buckling instability is still present. However, one still expects some effects of increasing attraction as E± is 
ramped up. One possibility is a structural transition to a different lattice geometry. In the remainder of the paper, 
we show that there is indeed a region of E±-n space where a square lattice becomes more favorable to the triangular 
lattice for the Wigner crystal. 

Bonsall and Maradudin's classic work 10 (referred to as BM from now on) has shown for a 2D electron lattice, where 
the electrons interact via a Coulomb force only, that the triangular (hexagonal) lattice is energetically the most stable 
of all five Bravais lattices. This is simple to understand physically, because for different 2D lattices corresponding 
to the same density, the triangular lattice is the one with largest lattice spacing; the electrons thus minimize the 
repulsive energy by staying as far away as possible from each other. Presumably, the triangular lattice is also the 
most stable for purely repulsive potentials of other forms, because the same argument holds. However, when one adds 
an attractive force, one of the other lattice shapes may become more favorable. For example, the lattice spacing is 
a factor of (2/v / 3) 1 ^ 2 smaller for a square lattice of the same density, and so the attractive energy can possibly be 
lowered by choosing this lattice geometry. Of course the number of nearest neighbors is also smaller for the square 
lattice, so which lattice geometry is energetically favorable depends on the exact forms and relative strengths of the 
attractive and repulsive interactions. 



IV. TRIANGULAR VS SQUARE LATTICE 

At finite temperature, the lattice geometry that is more favorable is the one with lower free energy F = E — TS. 
We will restrict ourselves to zero temperature, so that it is sufficient to consider the energy E of the two lattices. 

Using the two-electron potential (eq 4), one can simply sum up all pairwise interaction energies for the square lattice 
and the triangular lattice, and then compare. Instead, we will look at a slightly different quantity, as in BM 10 . We 
will consider the electron located at the origin and let E denote the energy of the interaction of this electron with all 
the other electrons. The total energy of the lattice of N electrons is then \NE. We compare triangular and square 
lattices: 

E™ = ]T T/(R™), E^ = £ V(R^) 

R,#0 R,#0 

Here R,™(SQ) runs over the positions of all the electron positions in the triangular (square) lattice, and V(r) is the 
potential (eq. 4) consisting of a Coulomb repulsion and a surface-mediated attraction. The two lattices each have the 
same density n. 

We will look at the difference, 

AE = £TR_£ SQ = AE coul + AE surf , (5) 
between the triangular and square lattices. The square lattice is more stable if AE is positive. The Coulomb part, 

(TR) e2 (SQ) e2 

Ai; cou i = T7TR - E ^SQ' ( 6 ) 

R,#0 1 R t #0 

is known to be negative, since the triangular lattice is the most stable under Coulomb forces alone, and has been 
studied in detail in BM 10 . In terms of density, their results are 

A^coui = (-3.921034)eV /2 - (-3.900265)eV /2 = -(0.020769)e 2 n 1/2 

As for the surface part, 



(SQ) 



AE sur( = e -f±{ - E K ° d R ™l AO + E ^° (|Rf Q | /lo) = e -f±S sur{ , (7) 

7r<T [ R;#0 R;#0 J 7r0 ' 

it is not a priori obvious that this is positive, but numerical calculations confirm that it is. The transition corresponds 
to the values of E± and n for which the two terms just cancel each other out, AE SUI f = —AE cou \, i.e., 



^transition = J (2^) (Q.020769) (2/ v^) 1 / 2 n 1 /2 (5 surf ) _1 . (8) 



Zero-Temperature Phase Diagram 
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FIG. 1: E± vs n phase diagram. Large dots, from numeric computation, map out the line where triangular and square lattices 
have the same energy. In the region above this curve, the square lattice is more favorable. The horizontal line indicates the 
spectrum instability, corresponding roughly to the transition to the surface-buckled, multi-electron dimple state. At higher 
densities this curve should move downward (eq 9); this density-dependence is not shown here. The inset shows the quantity 
Ssurf (eq 7), determined numerically, reaching a constant value at densities 3> Iq 2 . 



To map out the transition line exactly, we need to numerically calculate Ssurf; the numerical results are shown in 
figure 1. 

We can also analytically predict the dependence of S sur f and E± on the density for larger densities (n 3> Iq 2 ). The 
approximation K (r/l ) w — \n(r /Iq)6(1 — r) is good for r <C lo- Using this "logarithmic" approximation, we argue in 
appendix A that S SUI { becomes independent of n, at large n. Therefore from eq (8) one finds the shape Ej_ ~ ^/an 1 / 4 
for the transition curve in the E±-n plane, for densities much larger than Iq = pg/cr. 

Since the surface-mediated attraction falls off quickly for increasing inter-electron distances, one would need stronger 
electric fields at lower densities to have comparable attractive and repulsive interactions. Therefore the transition 
curve should rise steeply on the lowcr-n side, for n < Iq 2 . In the logarithmic approximation the left (lowcr-n) side of 
the transition curve is simply a vertical line at n = Iq . 



V. ZERO-TEMPERATURE PHASE DIAGRAM 



T he zer o-temperature phase diagram in the i?j_-density plane is shown in figure 1. E± and n are plotted in units 
of \J~g~/Tq and 1/Z§ respectively; this guarantees that the same phase diagram is applicable to different cryogenic 
substrates, one simply has to replace the numerical values of a and lo = \J cr/gp for the particular liquid being used. 

The transition line E±(n) is found by calculating S SUI f = — ^ TR -Ko(fi) + -^o( r i) numerically and then using 
eq (8). BM, in doing the corresponding sum for the Coulomb potential 10 , use an Ewald sum to convert the summation 
over 1/r to a summation over a faster-decaying function. Since K is rapidly-decaying itself, no such procedure is 
necessary for our case. Also, BM use a renormalization procedure (subtracting off contribution due to spread-out 
positive charge) to get a finite energy for each lattice geometry. This is again unnecessary because the sums converge 
for the short-ranged K$ function. 

Our analytical calculation of the high-density side of the phase-transition line is actually very good; a Ej_ ~ n 1 / 4 
fit matches the numeric curve quite spectacularly. On the lower-n side, the minimum of the curve is at about an 
order-of-magnitude lower than Iq 2 . The rise of the curve at small density is not quite as dramatic as the vertical line 
predicted by the logarithmic approximation. 
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At higher densities, there is an E± above which the electrons cluster into multi-electron dimples (which themselves 
form a triangular lattice) . As far as we know, a detailed calculation of the exact transition line between the single- 
electron Wigner crystal state and the surface-buckled multi-electron dimple lattice has never been performed, but a 
simple estimate is obtained by considering the instability in the spectrum of a uniformly charged liquid surface 8 : 

Ej_ = A^pga] 1 / 2 - {2ime) 2 . (9) 

The low-density part of this curve is the horizontal line in the phase diagram, figure 1. The rt-dependent deviation 
due to the (2nne) 2 term is substrate-dependent, even with our choice of units for E± and n, and is important only at 
very high densities, and so is omitted from the plot. 

VI. CONCLUDING REMARKS 

In summary, we have reported a novel effect of surface distortions on the Wigner crystal formed by electrons 
deposited on a liquid substrate. Our analysis shows that, at low electron densities, there is a significant portion 
of the zero-temperature E±-n phase diagram where a square lattice is energetically more favorable than the usual 
triangular lattice. Since we have not done a stability analysis of the square lattice under the combined interaction of 
eq (4), we cannot yet say whether the square lattice is stable, or whether some lattice geometry other than square 
and triangular is the actual stable geometry. However, the energy calculation proves quite clearly the presence of at 
least one structural phase transition at low densities. 

The structural transition in this system is particularly remarkable because it arises from the interplay of two 
very simple forces, a long-range repulsion and a short-range attraction. Other physical systems for which structural 
transitions have been discussed (crystalline solids 18 , flux-line lattices, skyrmion lattices 16,17 ) tend to involve far more 
complicated interactions between the constituents. 

Experiments on this system have, until now, probed only densities significantly larger than Iq~ 2 . For 4 Hc, the 
capillary length corresponds to l^ 2 w 400 cm~ 2 , while typical experimental Wigner-crystal densities arc in the range 
<~ 10 5 -10 9 cm~ 2 . The same situation holds for explored electron densities on other surfaces and interfaces. To the 
best of our knowledge, there have been no experimental efforts aimed at exploring the electron crystal structure at 
very low densities (n <~ Iq )■ One possible experimental signature of a structural transition to a different geometry 
would be a shift in the resonance positions in a Grimes- Adams 4 type experiment. 

The analysis in the present work, in addition to the prediction of at least one structural transition, poses several 
questions. First, there is the issue of the region of parameter space in which the triangular and square lattices are 
actually stable. Stability questions can be investigated by studying the dynamical matrix or the elastic constants 18 . 
Instability of both lattices in some part of the E±-n plane would indicate that a different geometry is more favorable 
than both the lattices we have considered. (Situations in which more than two lattice types are important have been 
encountered previously in the context of skyrmion lattices 17 .) In such a case there is the added question of what 
lattice, or other structure, is the actual stable one. One way to find out is to calculate the lattice energies for all 
five Bravais lattices (as done for the pure coulomb case in BM 10 ), or better yet, for all possible lattice geometries, 
parameterized in a suitable way. And finally, there is the question of crossover between single-electron lattice structures 
to many-particle dimple lattice structures as one increases the electron density at high E± . One immediately plausible 
speculation is that this transition may proceed via a dimerization. Several of these issues are the subject of ongoing 
calculations and will appear in a future publication. 
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APPENDIX A: DEPENDENCE OF S surf ON DENSITY 

For densities significantly larger than no — 1/£q, we use K (r/lo) « — ln(r/Z )#(?o — r ) an d get for the S sm t (defined 
in eq 7): 

(TR) (SQ) 

s smi « - Yl ln (l R ™lAo) + E ln ( 
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The summation for each lattice covers all the lattice points (electrons) within a circular area of radius lo, the number 
N = n(irl1) of electrons is the same for the two lattices. The two terms each contribute a term of magnitude ±-/Vln(Z ), 
which cancel. Rescaling (r-j = Riv}! 2 ), 

(TR) (SQ) 

S smf « - £ ln(r™) + £ ln(rf>) 

While each sum depends on the radius Zo" 1 ^ 2 of the circle, the difference does not. This can be seen by considering a 
radius r, and then increasing the radius by a small amount Sr; the number of electrons in the shell is SN = n x 2irr6r 
for either lattice, and therefore the change in each sum is equal to 5N logr, which cancel. There is thus no effect of 
increasing the radius Ion 1 / 2 of the circle we're summing over, i.e., S^urf is independent of the density n. A S sur f vs n 
plot is provided as an inset in figure 1. 
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